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NEUTRINO FLAVOR OSCILLATIONS IN ROTATING 

MATTER 

MAXIM DVORNIKOV tt 

Abstract. We study the evolution of the neutrinos system in rotating 
matter. Neutrinos are supposed to be mixed massive particles interact- 
ing with background fermions by means of the electroweak forces. First 
we find the solutions of wave equations for the neutrino mass eigenstates 
in matter. Then we study the behavior of neutrino flavor eigenstates in 
background matter. The problems of neutrino bound states and neu- 
trino flavor oscillations are discussed. We also derive the analog of the 
quantum mechanical evolution equation for the system of two flavor 
neutrinos in rotating matter and analyze its solution for the particular 
initial condition for neutrino flavor eigenstates. 
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1. Introduction 

VO 

' Nowadays it is acknowledged that neutrinos play a significant role in the 

£Sj ■ evolution of massive stars at the ultimate stages of their life. When the mass 

of a star is about 8 — 9 or 40 — 60 solar masses, such a star can explode as 
a supernova with the emission of great quantity of neutrinos carrying away 
almost 99% of the initial gravitational energy of a star [T]. 

Although neutrinos interact rather weakly with background matter, these 
particles are the key component in the dynamics of a supernova explosion. 
After the supernova explosion the core a massive star is converted into a 
compact dense object, a neutron star. Neutrinos are supposed to give the 
contribution to the subsequent evolution of a neutron star, i.e. causing its 
cooling [2 . Besides the direct influence to the supernova explosion process 
neutrinos can also affect various macroscopic characteristics of a neutron 
star. For example, great peculiar velocities of pulsars can be explained by 
the asymmetric neutrino emission J3j. It is also supposed that the emission 
of neutrinos can cause the spin-down of a rotating neutron star [4J . 
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The fact that neutrinos are massive particles has many important conse- 
quences. Unlike photons that immediately escape the region where they are 
created, neutrinos with relatively small initial energies can form bound states 
inside or in the vicinity of various astrophysical objects. The most proper 
candidates for the formation of such non-trivial states are relic neutrinos. 
The possibility of gravitational clustering of relic neutrinos was studied in 
Ref. [5]. It was discussed in Ref. [6] that both Dirac and Majorana neutrinos 
can also create a superffuid condensate due to the Higgs boson interactions. 
We considered the situation of gravitational trapping of neutrinos by a mas- 
sive black hole in Ref. [7J . The neutrino trapping in both curved space-time 
and rotating matter was studied in Ref. [8]. 

There is also a possibility that neutrinos emitted in a neutron star form 
bound orbits inside the star due to their collective interactions with neutron 
matter. This issue was studied in Ref. [9j- Neutrino trapping inside a 
rotating neutron star was discussed in Ref. |10j . This effect results from the 
neutrino electroweak interaction with inhomogeneously moving matter. 

Recently we developed an approach for the description of neutrino flavor 
and spin- flavor oscillations in various external fields The problem is 

formulated in terms of wave quantum mechanics and involves exact solutions 
to wave equations for a neutrino in an external field. In the present work we 
study the evolution of massive mixed neutrinos in inhomogeneously moving 
matter. In particular, we apply the treatment to the emission of neutrinos 
in a rotating neutron star. 

In Sec. [5] we give the general formulation of two neutrino flavors interacting 
with rotating matter of the type found in neutron stars. In Sees. [3] and 2] 
we find the solutions for the Dirac equation for a neutrino interacting with 
moving matter for massless and massive particles respectively. We also 
compare our solutions with the previously found ones \10\ I12j . In Sec. [5] 
we discuss the possibility for low energy neutrinos to form bound orbits 
inside a neutron star. Neutrino flavor oscillations in rotating matter are 
discussed in Sec.EJ In Appendix|A]we state the solution of the wave equation 
for a neutrino in vacuum in cylindrical coordinates. Finally in Sec. [7J we 
summarize our results. 



2. General formulation 

Let us first formulate the evolution of two neutrino flavor eigenstates, 
v\, A = a,/3, interacting with moving matter due to the exchange of the 
electroweak Z and bosons. Phenomenologically this interaction with 
matter can be implemented by means of a set of neutrino wave equations 
with the external fields /(* [13] shown below. In the flavor basis, neutri- 
nos also have a non-diagonal mass matrix (m\y). The Lagrangian for this 
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system is then: 

C= ^(i7^-/ A ^P L K- £ m AV^v, P L = ^(l-7 5 )- (1) 

\=a,/3 \\'=a,/3 

In our case of interest, the flavor a will be either [i or r, and the flavor /3 
will be e. 

The form of the currents are determined by the neutrino interactions 
with the medium. In the case of a neutron star the medium consists of elec- 
trons, protons, and neutrons, with number densities n e , n p and n n , respec- 
tively, and n e = n p , corresponding to electrically neutral matter. Therefore, 
for the standard model neutrino flavors a (either [i or r) and /3 = e, the 
corresponding external fields have the following expressions p2 

where Gf is the Fermi constant and 



/£ = -^#, fi = ^m-%)> ( 2 ) 



3e=( n e,n e v), j£ = (n n ,n n v), (3) 

are the hydro dynamical currents of each of the background fermion species. 
We also assume that all background fermions rotate as a rigid body, i.e., 
moving with same velocity v. 

To study the evolution of the system ([1]) we diagonalize the mass matrix 
(w-AA') by introducing the set of the neutrino mass eigenstates ip a , a = 1,2, 
with help of the matrix transformation: 

TT i (TT s [Uai U a2 \ (cose -sin6»\ 

V \ = V Uxatpa, (U\ a ) = = . . (4) 

a ^ 2 \U e i U e2 J \sm9 cosO J 

Here 9 is the vacuum mixing angle, where 9 = means v\ = v a and v 2 = v e - 
After diagonalization the Lagrangian (TjQ) reads 

£= ^Vfyt -m a )^ a - ^ dab^al^PL^b, (5) 

a=l,2 a,b=l,2 

where the matter interaction term contains the 2x2 matrix (g^ b ) in the mass 
eigenstate basis 

Gf /(2j£ sin 2 - # sin 20 \ 

l5ab j V sm 29 (2# cos 2 - j&) J ' W 

The Dirac equation for the neutrino mass eigenstates, obtained from 
Eq. , has then the form, 

(i 7 ^ - m a - Oft) 1>a - gfapPL A = 0, a^b, (7) 

where the last term corresponds to the off-diagonal elements of Eq. ([6]) , so it 
is an interaction with the medium that mixes the different mass eigenstates. 

We solve Eq. ([7]) treating the last term (the mixing of neutrino types) as a 
perturbation, so that at zeroth order the neutrino types are decoupled. Also, 
since the typical energy of a neutrino emitted in a neutron star is ~ lOMeV 
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whereas the neutrino masses do not exceed a few eV, the neutrinos are 
ultrarelativistic and we can treat the masses m a also as perturbations. 

Now, for a rigid rotating medium, the interaction depends on a velocity 
v(r) = fi x r, where r is the radius vector from the star center and fl is the 
angular velocity of the star. Therefore we define the positive potentials as: 

v __o _ \G F (n n - 2n e sin 2 0)/y/2, for a = 1, 
9aa ~ [G F (n n - 2n e cos 2 9)/y/2, for a = 2. 

Note that V a > in Eq. ([8]) since n n > n e in a neutron star and g® a < 
(a = 1,2). 

In order to proceed with the description of the evolution of the system ([7]) 
we should have the energy levels and wave fuctions of the mass eigenstates 
tp a which correspond to the wave equation ([7J) at the absence of the mixing 
term g^ b . These quantities will be found in Sees. ElandHI 



3. Solution of the wave equation for a neutrino in rotating 
matter in the limit of zero mass 

In this section we find the solution of Eq. ([7]) at the absence of the mixing 
between different mass eigenstates due to the interaction with matter, i.e. 
we put the coefficient g^ 2 to zero. 

This case corresponds to a single unmixed neutrino interacting with an 
external axial-vector field. The wave equation ([7]) is transformed to the 
form PUQ2], 

(i 7 ^ - m - g^P L )^ = 0. (9) 

Here we are interested in the case of ultrarelativistic neutrinos, i.e when the 
mass in Eq. is much smaller than the energy. 

The equations of motion for the left-handed rj and right-handed £ chiral 
components of the spinor t/j T = (£, rj) decouple in the m = limit. The mass 
contribution can be included in perturbation theory (see Seed]). Using the 
chiral basis for the 7^ matrices in the convention of Ref . [16] , 

MW). 

the Dirac equation for the left-handed component r\ of the neutrino is 

if) = {ia ■ V + (T^}??, (11) 

where the matter interaction term is cr^g^ = go + cr ■ g. 

The external field g^ with non-zero spatial components corresponds to the 
external electroweak field due to the moving background matter. In analogy 
to Eq. ([S]) we obtain for the three-vector part, g = gov. We want to study 
the neutrino states inside a neutron star rotating with the angular velocity 
fi directed along the z axis, ft = Qe z . As we mentioned in Sec. [21 since 



MAXIM DVORNIKOV : NEUTRINO FLAVOR OSCILLATIONS ... 5 

go < in the case of a neutron star, we redefine this potential as V = —go- 
Accordingly, 

g = Vn(ye x - xe y ). (12) 
It is natural to use cylindrical coordinates (r, <ft, z) to solve Eq. (jlip with 



kr • V =i 



8 Z e-^ldr - (i/r)fy] 

j*[d r + (i/r)fy] -8 Z 



Taking into account that Eqs. (jlip and (|13p do not depend on z, we look for 
a stationary solution of the form rj ~ u(r,cj))e~^ Et ~ PzZ \ Using the analysis 
of Ref. [T7] we can further separate the two coordinates r and <j) in the 
two-component spinor u using auxiliary functions F\ 2 (p) as 



where p = VQr 2 , and / is an integer so that the function of (f> is single- 
valued. Since the system is invariant under rotations around the Z axis, u 
has to definite value of J z (total angular momentum), which in our notation 
is equal to / — 1/2. 

The functions F\ t2 satisfy the coupled equations: 

i^JVOTp (2d p - 1 - F 1 = i/fcFi = {E-p z + V)F 2 , 

i^/Vnp~ (2d p + 1 + -\ F 2 = iR 2 F 2 = (E + p z + V)F 1 , (15) 

which can be separated into two independent second order differential equa- 
tions: 

d 2 d I p (l-l) 2 ) 

* + f + „_Lil_£ pj^.o. (it) 

dp 2 dp 2 4 4pJ 

where the parameter k is related to the energy in the form, 

Let us just solve Eq. (fT6|) for the function F\. Eq. (fTTj) can be solved for 
F2 by analogy. Expressing Fi(p) = e~ p / 2 p( l ~ 1 ^ 2 u(p), the new function w(p) 
obeys an associated Laguerre equation, 

pu" + (/ - + (k - l)u = 0, (19) 

whose solutions are the associated Laguerre polynomials u{p) ~ Q l s r 1 (p). 
Here s = k — I is the radial quantum number, and k, defined in Eq. f)18[) . is 
related to the neutrino energy. The function F\ is then a Laguerre function, 
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F\(p) = The Laguerre functions In,s(p) & r e denned in terms of 

the associated Laguerre polynomials Q l s (p) (where n = s + I) as 

I s+ljS (p) =^ = L=e-''/y/ 2 Qi(p), 

Q l s {p) =ePp- l ^- s {p s+l e'P)- (20) 

Another common definition of the associated Laguerre polynomials is L l s (p) = 
Q l s {p)/s\ The Laguerre polynomials satisfy the recursive relation, 

Q l s - l {p) = Q l s {p)-sQ l s - l { P ). (21) 

We also mention that I n ,s(p) an d Q l s {p) satisfy the integral relations, 

;>oo 

/ I^ s (p)I n -i. s (p)^/p dp = Vn, 
Jo 

POO 

/ e~Pp l Q l s (p)Q l Ap) dp = 5 ss ,s!(l + s )! (22) 
Jo 

In order to have normalizable functions at the origin, I must be a non- 
negative integer. Also, the solution diverges at large radii, unless s is a 
non-negative integer. Therefore, in order to have well-behaved solutions, 
the smaller values of s and k must be integers, in which case the lower 
values of energy, E = — V ± \J AVQk + pi, are discrete: 

k — y n, E -> E n = -V ± y/Wn n + p 2 z . (23) 

On the other hand, for large enough energies the divergent behavior of the 
solution falls outside the star radius, in which case k has no restrictions and 
it becomes a continuous variable. 

Notice that the energy levels in Eq. ([23]) are different from the analogous 
expressions obtained in Ref. |10| . E = — V ± y / 2FJ7n + pf. It was claimed 
in Ref. [10j that there is an analogy between the charged particle dynamics 
in an electromagnetic field and the neutrino motion in matter. However 
this analogy is just superficial. An electromagnetic field is gauge invariant. 
Therefore when we study the motion of an electron in an external magnetic 
field B = (0,0, B), we can choose any gauge. In a gauge that keeps the 
cylindrical symmetry explicit, the vector potential A = (—yB/2,xB/2,0) is 
required. The analog of this gauge is adopted in the present work. If the 
electron motion is studied in cartesian coordinates, the Landau gauge is more 
convenient, A = (0,xB,0). This kind of gauge was used in Ref. [TO]. As 
shown in Ref. |17j . in the electromagnetic problem both gauges must give the 
same energy spectrum for the electron. However for the motion of a neutrino 
in rotating matter, the situation is different: there is no gauge freedom in 
this case. If one uses the analog of the Landau gauge, as in Ref. [ID], one 
underestimates the matter contribution to the dispersion relation. 
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Finally, to derive the lower component of the neutrino spinor u, we just 
need to use the identities (see Ref. [T7]). 

RlIn-lM = -VWn^I n , s ( P ), R 2 InAp) = y/Wftoln-iM, (24) 

to find the function i^O )- Thus we get the complete two-component spinor 
in the form 

u(r 6) - ^ (CJa-iM^-W} m 

u{rA) ~ V ic 2 i n , s ( P )^ )■ (25) 



The coefficients C\^, are related to each other due to Eq. (fT5|) as 

VAVttn Ci + (E + V - p z )C 2 = 0, (26) 

and their norm can be chosen to satisfy Cf + Cf = 1. 

Let us now discuss the limit of small angular velocities to establish the 
connection with the non-rotating case. For simplicity we study the case 
p z = which is examined in Sec. [5J The limit f2 — > should be taken 
together with n — > oo, so that $7 • n = constant. Using the identity |18j . 

lim I nin _,([p ± r] 2 /4n) = J,( P± r), (27) 

n— >oo 

we reproduce the neutrino wave functions in vacuum (|56p . Therefore one 



can identify V 4US7n with the neutrino momentum in the equatorial plane, 
p±, in Eq. ®. 



4. Approximate solution of the wave equation for a massive 
neutrino in rotating matter 

In this section we will study the effect of the rotation of matter on the 
single massive neutrino without mixing, analogously to the treatment of 
Sec. El However, now we will take into account the contribution of the 
neutrino mass to the energy levels (|23p using the perturbation theory. 

For a massive particle one should take into account both spinors £ and rj 
in Eq. ([9|). The coupled wave equations for these spinors have the following 
form: 

i£ =(er • p)£ - mry, 

vfj = - (a ■ p)rj - m£ + [g + (a ■ g)]r], (28) 
where the vector g is defined in Eq. (|12p . 



Looking for the stationary solutions of Eq. (|28p . £ ~ e~ lEt and r\ ~ e~ lEt , 
and excluding spinor £ from Eq. (|28p we get the only differential equation 
for the spinor r\ , 

[E 2 - m 2 - p 2 + VE - E(o- ■ g) - V(a ■ p) + (a ■ p)(cr • g)]^ = 0. (29) 



Note that one should take into account the non-commutativity of the oper- 
ator p and vector g since the latter depends on the spatial coordinates. 
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As in Sec. [3] we will use cylindrical coordinates (r, 0, z) to analyze Eq. ([29]) . 
It is convenient to rewrite this equation for each of the components of the 
spinor rj T = (771,772), 

1 d 2 

E 2 -m 2 + VE + 8 2 + -d r + ^ + d 2 + iVVtd^ 



+ iVd z + Vn(rd r + 2) } m 



-e-^V < i 



n r (iE-d z )} m 



d r - -84. 



\ d 2 

E 2 -m 2 + ve + d 2 + -d r + ^ + d 2 + ivnd 4 



ivd z -vn(rd r + 2)} m 



&V < i 



+ Qr(iE + d z ) \ 771 



d r + -d<k 
r 



We look for the solution of Eq. (|3U|) in the following form: 

1 



77 = e rpzZ u(r, (f)), u(r, 



e i(/-l)<£ p/ r > 

W V e il ^F 2 (r) 



(30) 



(31) 



where -^1,2(7") are the new unknown functions [see Eq. (|14|) ]. From Eq. (|3U|) 
we derive the equations for the functions Fi t 2(r), 



E 2 -m 2 -p 2 + VE + d 2 + -d r 

r 

- Vp z + Vn(rd r + 2)\F 1 

r 



-iV 



d r + 



nr(E- Pz )}F 2 , 



1 l 2 

E 2 -m 2 - p 2 z + VE + d 2 + -d r - - 



VQl 



+ V Pz -Vn(rd r + 2) }F 2 

(I -I) 



-iV 



dr 



+ nr(E + p z )}F 1 . 



(32) 



Introducing the dimensionless variable p = Vilr 2 , as in Sec. [31 and using 
the properties of the operators R12, defined in Eq. (fT5j) . 



i?ii? 2 =wn P d p + a. 



r 2 Ri =4vn ( pd p + d p 



1 -1 P i 2 

1 P (i -if 



2 4 



(33) 
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we rewrite Eq. (|32p in the following form: 



{E 2 -m 2 -p\ + VE- Vp z + 3Vn + R 2 Ri + y / VQpR 2 }F 1 



= -i [VR 2 - ^/VVTp{V + E - p 2 )} F 2 , 



{E 2 -m 2 - P 2 Z + VE- Vp z - 3Vn + R X R 2 - ^fVVLpR 1 }F 2 



= -i [VR X + y/VQp'(y + E + Pz )}f 1 . 



(34) 



To study the contribution of the neutrino mass to the neutrino energy 
spectrum (|23|) we discuss the situation of the neutrino bound states and 
take into account neutrino mass with help of the perturbation theory. It 
means that in Eq. (|34[) we can use the expressions for the wave functions F\ t2 
presented in Eq. (|25p which correspond to the massless neutrino: Fi(p) = 
Ci/ n _i iS (p) and F 2 (p) = \C 2 I niS (p). For this kind of wave functions we get 
from Eq. (|34[) the following relations: 



{Ct[E 2 -m 2 -p 2 z + V(E-p z + 40 - 40n)] - C 2 VVwnn}l n - hs (p) 
+ ^/VQp~(V + E - p z )C 2 In,s(p) + 2Vn y /p{n-l)C 1 I n - 2 ,s{p) = 0, (35) 

{dVVwnn - C 2 [E 2 - m 2 - p\ + V(E +p z - 40 - 4fin)] }l n , s (p) 
- ^/vfTp{V + E + pJdln-iM - 2VQ^/p(n + l)C 2 I Q+1 , s (p) = 0. 



To obtain Eq. ()35p we use the known properties of the Laguerre functions, 



which can be found in Ref. [19j . 

Let us study the situation when E 2 3> VQ. and n 3> 1. The former 
condition is always satisfied for any realistic situations. Indeed, for a neutron 
star with n n = 10 38 cm -3 we get that V ~ 10 eV. If we suppose that the 
energy has the value of several electron- Volts (we will see in Sec. [5] that a 
bound state can be formed for such low energy neutrinos) and a neutron star 
has the angular velocity f2 = 10 3 s _1 (~ 10~ 13 eV), we get that the condition 
E 2 3> V£l is satisfied. The latter condition, n 3> 1, is also valid (see Sec. [5]). 
It is received there that the critical value of the quantum number n at which 
a bound state is still possible is equal to 10 11 . It means that this condition 
is fulfilled. 

Taking into account the approximations made above we obtain that the 
energy spectrum has the form, 



However, instead of the relation (|26p . which is valid for the massless case, 
we have the corrected one, 




(36) 



E = -V± s/Wttn + pJTm?. 



(37) 



VWQn C x + (E + V + 40 - p z )C 2 = 0. 



(38) 
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In Sec. Owe obtained that the quantity p e g = \J 4Vfin + p 2 z has the meaning 
of the effective momentum of a neutrino. Therefore we get in Eq. (j37[) that in 
the limit n 3> 1 the neutrino energy receives the correction m 2 /2p e s due to 
the non-zero mass, as it should be. Note that this energy spectrum coincides 
with the analogous relation found in Ref. [12J. 

5. LOW ENERGY WITHOUT FLAVOR MIXING: BOUND STATES 

In Sees. [3] and H] we found the basis spinors ([25]) and the energy levels ([571) 
of a single neutrino mass eigenstate interacting with rotating background 
matter. Now we apply these results for the description of the evolution of 
the system ([7]). 

Given the axial symmetry of the problem, we use cylindrical coordinates 
and the eigenfunctions = (Ca^Va) of Eq. ([7J must be of the form r\ a ~ 
rj a (r, <fi)e~*(® at ~ PzZ ' . We find the energy eigenvalues in terms of an index "n" 
[see Eq. (JSTJ], 

E {a)± = _y a ± v /4y a n n + p 2 +m 2 ) n = , 1, 2, . . . (39) 

Here E^ + is the energy of a particle (neutrino), containing an attractive 
potential "— V a " , while the negative value E^ must be understood as 
—E^ , the positive energy of an antiparticle (antineutrino), containg a 
potential term a +V a " , which is repulsive. We included the neutrino mass 
m a in the above expression for the energy, although we will neglect it in the 
spinors. 

In a rotating medium, the effect of rotation is largest for neutrinos prop- 
agating in the equatorial plane. Choosing Z as the rotation axis, we then 
look for solutions which are z-independent, i.e. p z = 0. 

The corresponding two-component spinors are given in terms of Laguerre 
functions I njS (/j) with an energy index "n" and a radial index "s" [see 
Eqs. (|25D and (1261) ]: 

^ ( ^ ) = klTil4^ J' Z = n " S ' (40) 

and where p a = V a VLr 2 , a = 1,2, is a dimensionless radial coordinate. For 
further details of the derivation of Eqs. ([39p and ([40 p the reader is referred 
to Sees. [3jand[U 

In this case, when "n" is a non-negative integer, the energies in Eq. ()39p 
are discrete values and the wavefunctions decay to zero as r — > oo. This is 
a required condition for low enough energies, see discussion below Eq. (j44p . 
otherwise the wavefunctions would diverge inside the star before reaching 
the star radius. 

Instead, for neutrinos with larger energies, the wavefunction reaches the 
edge of the star and should continue outside. For those cases the neutrino 
energy and the index "n" become continuous variables (we change the name 
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n — } k in this case). The solution inside the star now has the more general 
form: 



(41) 



( c[ in) e-^ 2 pt m F(l - k, I, Pa^-W/i (I - l)!v^A 

where F(a, {3, z) is a confluent hypergeometric function and the coefficients 
satisfy the relation 

y/WatlK C? n) + (E a + V a - Pz ) Ct ] = 0, (42) 

considering the energy E a with the same expression as in Eq. (|39p , but with 
n— > k, a continuous variable. We derive Eq. (|42p from Eq. ()26[) changing 
the norm of the coefficients C\ 2 there as: 



cS^C^Jj-^—, (43) 



(n-0 



and using the property of the confluent hypergeometric function, F(l — n, 1 + 
l,p) = l\Q l n _ l ( P )/n\. 

For the wave function outside the star one must use the outgoing wave 
solution in vacuum given in Eq. (|63p . i.e. 

where -ff/ 1 ^ are Hankel functions of the first kind and p± = \J E 2 — — m 2 
is the momentum perpendicular to the rotation axis. 

The coefficients cj 1 ^ and c[°2^ are related due to the continuity of the 
wave functions ([IT]) and (|4"4"|) at the neutron star surface, Ua^(R,(j)) = 
u^ ut \R,4>). Eqs. (|4"2"jh ([59)) and ([66]) completely define the coefficients for 
the solution corresponding to an unbound wave function. 

As mentioned above, neutrinos could form bound states inside a neutron 
star, provided the energy is small enough. In those cases the energy (j39j) 
assumes discrete values, because otherwise the wave function would diverge 
before reaching the edge of the star. In Fig. [JJ we present an example of 
wave function [see Eq. ()40p ] for I = 10 and s = 15 (solid line). It can be seen 
that the solution is a rapidly oscillating function in the inner regions and 
approaches zero for large radii. We also present in Fig. [JJ the wave functions 
for I = 10 and s = 15.1 (dashed line) as well as for I = 10 and s = 14.9 
(dash-dot line), showing the divergence at large radii. 

From the fact that can be negative for a medium dominated by 



neutrons, as seen in Eq. (|39|) . neutrinos of low kinetic energy can form 
bound states inside the star. On the other hand, antineutrinos cannot be 
bound because in their case the potential is repulsive. Limiting the analysis 
to p z = only, the energy spectrum for the bound states is discrete. The 
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X 

Figure 1 . The solutions of Eq. (JTUJ) for I = 10 and s = 15 
(solid line), I = 10 and s = 15.1 (dashed line) as well as 
I = 10 and s = 14.9 (dash-dot line). The function I n , s (x) 
at non-integer "s" should be understood as the confluent hy- 
pergeometric function F(k — 1, 1 + 1, x) in Eq. ()4ip . 



maximum value of "n" for bound states can be estimated with the equation 

E { n> + = 0, 

(45) 

This condition is independent from (but consistent with) the condition that 
the wave function should not diverge before reaching the star radius. The 
position of the last maximum in the wave function is approximately at p CT it ~ 
2(n + s), and this point should be inside the star, i.e. p cr i t < V a QR 2 , where 
R is the star radius (remember the definition p = V£lr 2 ). This condition 
means n < V a ftR 2 /2. 

We can estimate the index "n" in these bounds for the case of a neutron 
star. Taking a reference value for the neutron density n n = 10 38 cm -3 (close 
to nuclear density), we get a potential V a ~ 7 eV. This is much larger than 
the natural neutrino masses expected from oscillation experiments, so the 
mass in Eq. (j39|) and Eq. ()45|) is indeed negligible. Now, for the angular 
speed Vt, W6 CcLn take the Keplericin angular velocity ^max — 

a/ GM/R 3 (see 

Rcf. |20j) as a rough upper limit. Taking M = 1.5M Q and R = 10km one 
obtains for f2 max ~ 10 4 s" 1 (typical values are two or more orders of magni- 
tude smaller). Eq. (|45p then gives n ~ 10 11 , a rather large number of bound 
states as well as nodes within the star. On the other hand, the condition 
for the wave function to fade before reaching the edge, n < V a £lR 2 /2, gives 
also n ~ 10 11 . 

The particles inside the neutron star can be created in the form of neutrino- 
antineutrino pairs [2]. The fact that low energy neutrinos can be trapped 
by the star and antineutrinos will be accelerated out of the star makes it 
possible the existence of the low energy antineutrinos luminosity. 
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The validity of the wave function method has one main drawback con- 
cerning wave coherence: it must be assumed that incoherent scattering can 
be neglected all over the star. This is particularly unlikely if the number 
of nodes is as large as 10 11 . Another drawback of this treatment concerns 
neutrino production: most neutrinos are produced with energies much larger 
than a few electron- Volts, and again, the production is in general an inco- 
herent process. 



6. Effect of rotation on flavor mixing of low energy 

neutrinos 

So far we have neglected the off-diagonal elements of Eq. ([6]) which mix 
the neutrino mass states. Let us now include them as perturbations to study 
the effect of rotation in the mixing of neutrino flavors. The general solution 
of Eq. ([7]) can be stated in the following form (see also Ref . [11] ) : 

oo 

Va (r, 4>, t) = (ail <ns(r, <t>) exp[-iE^ + t] (46) 

n,s=0 

+ b<£ u- ns (r, <f>) exp[-L4 a H]) , a = 1, 2. 

In Eq. (|46p we consider the full state as an expansion in the eigenstates 
with time-dependent coefficients tins (t) and b$ it) , requiring the expansion 
to satisfy the evolution equation Q, and imposing the initial condition that 
the neutrino wavefunctions at t = are of flavor f3 = e only. 

Substituting the solution into Eq. and using orthonormality of the 
spinors in Eq. (|40p we obtain the following set of ordinary differential equa- 
tions for the functions a^J (t) : 

4/nHt)= £ {/<L(r)(^ 6 ^)<nvWd 2 r} 

n',s'=0 ^ J J 
+ lj <L(r)(5ab^)%nv( r ) d2r } 

x exp[i(4 a )+ - E®~)t] b%(t), a / b, (47) 

and a similar set of equations for bnJ(t). 

Now, since we are interested in the neutron star matter, we can assume 
densities close to nuclear matter or above, dominated by neutrons (n e <C n n ). 
Within this approximation, the spinor products in Eq. (|47p can be easily 
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calculated, 

J <L(r)(< b ^) u M' s '( r ) d2r = /\L( r )(i^)v S 'W d2r 

~ 5i2 <^Z'<W —5i V (25 ss iy/n- Vs + l<5 s , s '_i - \/s<^,s'+i)> 

y <L( r )(^^) u£ nV (r)d 2 r = /<ns(r)(^^)< nV (r)d 2 r 

« -^y^Mv^V+i - >ATTV-i)> ( 48 ) 

where Z = n — s, Z' = n' — s', g® 2 is the potential that mixes neutrino 
types, given in Eq. ©, and V = GFn n j\pl is the approximation of V\ ~ Vi 
neglecting the electron density n e . Here we must keep terms linear in n e , 
otherwise g® 2 would vanish, and the case of mixing would become trivial. 
To linear order in n e , we can take the arguments of the Laguerre functions 
in Eq. ((ID]) to be equal, p\ = P2, neglecting n e in Eq. ([8]). To derive these 
results we used the properties ([2T]) and ([22]) of the Laguerre functions, given 
in Eq. (p2]) . 

The time evolution for a£j(t) given in Eq. (|47p . which is due to the fla- 
vor off-diagonal perturbation, has contributions from many levels {n, s}, 
making a general analysis rather complicated. One simplification occurs if 
we disregard neutrino-antineutrino creation or annihilation, which means 
no transitions between u + and u~ states. In what follows we disregard 
antineutrinos, so we neglect the coefficients bus (t) altogether. Another sim- 
plification arises if we consider wave packets very narrow in energy, so that 
only the nearest states are involved. A final simplification occurs in the case 
s > I, where level transitions are negligible, thus reducing the problem to 
solving a two-state quantum system. 

The situation s 3> Z, occurs when the neutrino angular momentum at 
emission is small, corresponding to wave functions which are large near the 
center of the neutron star. In this case the coefficients a$ are, in practice, 
functions of just one quantum number s, since s = n — Z while Z is negligible. 
Assuming depend smoothly on s, i.e. a{ a ^ « &i±i; and usin § Eq. ([39 
we get an evolution equation of the form, 




A 

oj/2 




afe 1 ^ 2 , (49) 



where 



u =E^ 1)+ - Ef )+ = V 2 -Vi + yJW^ls + m\- J4V 2 ns + mf , 

A =g° 12 . (50) 

Since, for non-bound states, s is a very large and continuous variable, we 
can define a continuous "momentum" variable inside the medium, p e g = 
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y/4V£ls, chosen to be the same for both neutrino types, as it is usually done 
in the treatment of non-rotating media |21| . In terms of more conventional 
quantities, A and uj are then 

uj = V2GYn e cos 29, A = —=n e sm29, (51) 

2f> c fT V2 

In Eqs. (|50|) and (|5ip we take into account the non-zero masses of the neu- 
trino mass eigenstates mi t 2 and use the standard notation 5m? = — m|. 
As one can see, Eqs. (|i9|) and ([51]) are practically independent of the rota- 
tion velocity Q, and so the evolution equation reduces to the known case of 
neutrino oscillations in non- moving background matter. 

For the rotation to have any significant effect on neutrino flavor oscilla- 
tions the linear velocities of the matter motion should reach very high values. 
For the realistic angular velocities of pulsars ~ 10 3 s" 1 and a neutron star 
with radius of ~ 10 — 20 km the matter velocity can be about 0.1, which is 
not enough to significantly change the transition probability. Therefore the 
effect potentially interesting from the phenomenological point of view is the 
trapping of low energy neutrinos discussed in Sec. [5j 



7. Conclusion 

Summarizing we mention that we studied neutrino flavor oscillations in 
rotating matter. The analysis was carried out in frames of the relativistic 
quantum mechanics. We used the exact solutions of the wave equation for 
a neutrino weakly interacting with inhomogeniously moving matter to solve 
the initial condition problem for the system of mixed flavor neutrinos. Note 
that the use of the relativistic quantum mechanics method |11| is essential 
in describing neutrino flavor oscillations since it takes into account the co- 
ordinate dependence of the neutrino wave functions (flD"|) . It was possible to 
derive the Schrodinger like evolution equation for the two component neu- 
trino wave function for the important case of neutrinos with small angular 
momentum. This situation corresponds to the particles emitted close to the 
central region of a neutron star. It was found that rotation does not change 
the dynamics of neutrino flavor oscillations, i.e. the quantum mechanical de- 
scription of neutrino oscillations is insensitive to the rotation of background 
matter. 

We have studied the possibility for the existence of neutrino bound states 
inside a neutron star. The cases of neutrinos and antineutrinos are differ- 
ent. We revealed that low energy neutrinos can be trapped by the rotating 
neutron star whereas antineutrinos always escape. The applicability of the 
relativistic quantum mechanics method was analyzed. Despite this approach 
allows one to account for the interaction with external fields exactly it can- 
not be used for the description of the evolution of high energy neutrinos in 
dense rotating matter since these particles experience many incoherent colli- 
sions and it is difficult to form the coherent cylindrical wave inside a star. It 



16 



APPL. COMPUT. MATH., VOL. 9, NO..., 2010 



means that the relativistic quantum mechanics is applicable for later stages 
of the neutron star evolution when one can neglect the multiple neutrino 
collisions with background matter. 
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Appendix A. Solution to the wave equation for a neutrino in 

VACUUM IN CYLINDRICAL COORDINATES 

In this Appendix we find the solution of the wave equation for a neutrino 
in vacuum in cylindrical coordinates. In the chiral basis for the Dirac 7- 
matrices [see Eq. (jlOp ] the Dirac field i^(x) is expressed in terms of two 
2-component chiral fields, rj(x) and as ^ T = (£, 77). In cylindrical 
coordinates (r, 4>, z), the spinors £ and rj satisfy the differential equations 

ii =(tr • p)£-7W7, 

if] = — (<r • p)rj — m£. (52) 



We look for stationary solutions of Eq. (|52|) . which must be of the form, 

£(r, t) = e-^ Et - p ^w(r, 0), r](r, t) = e-^ Et - p ^u(r, 0), (53) 

where E is the energy and p z is the z component of the neutrino momentum. 
We can separate the variables (r, eft) introducing radial functions Gi %(r) and 
Fi,2(r), 

Wl (r, 4>) =^=e i V- 1 ^G 1 (r), w 2 (r, cp) = -±=j l *G 2 {r), 

V Z7T V 27T 



u l {rA)=^=e^- 1 ^F l {r), u 2 (r,0) = -^=e w *F 2 (r) ) (54) 

V Z7T V 27T 



where Z measures the z-component of the orbital angular momentum. For 
the radial functions G\ t 2 an d ^1,2 we obtain the system of the differential 
equations, 

(E - Pz )Gx +i(d r + l -}G 2 = - mF u 
(E + Pz )G 2 + i(d r - —) Gi=- mF 2 , 



{E+pjFx -i[d r + -)F 2 = - mGx, 



(E-p z )F 2 -i[d r -—\F x =-mG 2 . (55) 
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The solutions of Eq. (|55p can be given in terms of Bessel functions. For 
standing waves that are regular at the origin, the functions w(r, 4>) and 
u(r, 4>) are 

. 1 fBi < 7,_ 1 (pj_r)e i ( , - 1 >* 

w{r,<p) — ' 



<r,<P) = ^r:^!lu, ), (56) 



V iB 2 J l {p L r)e ll <t> 

Yf V iC 2 Ji(p±r)e 

where = -E 2 — p 2 — m2 ) -^1.2 an d Ci,2 are the undefined coefficients. 
Using the identity for the Bessel functions, 

J'l(x) + -Ji{x) = Ji-i(x), J , l _ 1 (x)- l —^Ji- 1 (x) = -Ji(x), (57) 

and Eq. ([53]) we find that the coefficients .6^2 and obey the system of 
the following algebraic equations: 

(E-p z )B 1 -p ± B 2 =-mC l , 
{E + p z )B 2 - p L B x = - mC 2 , 
(E +p z )C 1 +p ± C 2 =-mB 1 , 

(E-p z )C 2 +p L C 1 = -mB 2 . (58) 

In the limit of the small neutrino mass one can see that the equations for 
Bi t 2 and C\ :2 decouple giving one relation only for the coefficients C\ j2 , 



y/E+ Pz d + tJE - Pz C 2 = 0. (59) 
Normalizing the wave function as 

d 2 r u E ,i(r, ^ u E ', v (r, 4>) = E6(E - E')6 W , (60) 



and using the integral property of Bessel functions 

1 

drrJi(kr)Ji(k'r) = -5(k - k ), (61) 

one finds that the coefficients C\ % also must satisfy the relation 

Cl + Cl = E 2 . (62) 

Instead of standing waves like Eq. (|56|) . we can have outgoing waves. 
These solutions are similar, but in terms of Hankel functions of the first 
kind, h\ 1] {x), 

J_ (CxH^p^e^ 1 - 1 ^ 
f^ni \ \C 2 H[ l \p L r)e li 

The Hankel functions have the asymptotic behaviour of outgoing radial 
waves in the two dimensional plane 

T 



h\ 1 \x) « \j — exp 
" irx 



( irl 









(64) 
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Instead of the normalization (|60p one has to define the particle flux at large 
distances, j^. For a radial flow in two dimensions, the flux of particles 
through a circle of the large radius r has the form, 

joo = lim / rdcj) vJ (r, (j>) u(r, (f>), (65) 

r-J-oo J Q 

where we have assumed relativistic particles. From this expression, one can 
derive norm for the coefficients C± t 2, 

C! + Ci = ^^. (66) 



On the other hand, the relation (|59|) is still valid for the wave function (|63|) . 
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